Damping of Collective Nuclear Motion and 
Thermodynamic Properties of Nuclei beyond Mean 

Field** 

Hong-Gang Luo 1 ' 2 , W. Cassing 1 and Shun- Jin Wang 1,2 

1. Institut fur Theoretische Physik, Universitat Giessen 
35392 Giessen, Germany 

2. Department of Modern Physics, Lanzhou University, 
m ; Lanzhou 730000, PR China 



February 9, 2008 

> 
m 

Abstract 

The dynamical description of correlated nuclear motion is based on a set of 
, coupled equations of motion for the one-body density matrix p(ll';t) and the 

two-body correlation function C2(12, l'2';t), which is obtained from the density- 
matrix hierarchy beyond conventional mean-field approaches by truncating 3- 
body correlations. The resulting equations nonperturbatively describe particle- 
particle collisions (short-range correlations) as well as particle-hole interactions 
(long-range correlations). Within a basis of time-dependent Hartree-Fock states 
these equations of motion are solved for collective vibrations of 40 Ca at several 
finite thermal excitation energies corresponding to temperatures T = — 6 MeV. 
Transport coefficients for friction and diffusion are extracted from the explicit 
solutions in comparison to the solutions of the associated TDHF, VUU, Vlasov 
or damped quantum oscillator equations of motion. We find that the actual 
magnitude of the transport coefficients is strongly influenced by partlicle-hole 
correlations at low temperature which generate large fluctuations in the nuclear 
shape degrees of freedom. Thermodynamically, the specific heat and the entropy 
of the system as a function of temperature does not differ much from the mean- 
field limit except for a bump in the specific heat around T ~ 4 MeV which we 
attribute to the melting of shell effects in the correlated system. 
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1 Introduction 



The damping of nuclear collective motion and the role of short-range and long-range 
nuclear correlations Jl|, || in the relaxation still is left as an open problem from the 
fully microscopic point of view. Whereas at intermediate energy heavy-ion collisions 
the system is essentially damped by in-medium nucleon-nucleon collisions within the 
time-dependent mean field JJ], ^, [IJ , the relaxation of giant resonances at low excitation 
energy appears to be dominated by residual 2 particle - 2 hole (2p-2h) correlations || [7|] 
that are not incorporated in standard transport theories || f|, Though the latter 
approaches have proven to quite reliably describe one-body observables in heavy-ion re- 
actions, they fail in describing fluctuation phenomena like 'multifragmentation'. Thus 
Langevin forces - that are related to fluctuations in the particle collision number or in 
momentum space - have been proposed in addition to the Vlasov-Uehling-Uhlenbeck 
(VUU) dynamics || - [JI8] . Since the various approaches involve quite different as- 
sumptions about the dynamical origin and the size of 'fluctuations', a more rigorous 
analysis within the framework of a nonperturbative quantum two-body theory appears 
necessary. 

More specifically, nuclear giant resonances of hot nuclei have been studied in detail 



in both experimental (19], [21], ^2j and theoretical nuclear physics |23|] - j33] to learn 
about the nuclear collectivity and the lifetime of these modes. The stability of the 
damping widths against temperature, the effects of shape degrees of freedom on the 
damping process, the role played by particle-hole (long-range) correlations, and the 
distinction among the damping processes in different dynamical models have presented 
open and interesting problems. Though several authors have addressed these problems 
within various approaches [|, % 0, |2|, |25], g§, [27], [28], [29], |32|, [E| 0, g§ © M , many 
microscopic aspects remain unclear. It is the purpose of this paper to investigate these 
problems from a systematical point of view and in a more detailed way especially in a 
comparison of quantum mechanical and semiclassical approaches. 

In Section 2 we briefly review the underlying microscopic approach obtained from 
the correlation dynamical truncation scheme ]38|] and present the approximations per- 
formed for the actual calculations. Section 3 is devoted to a presentation of the micro- 
scopic results in case of isoscalar quadrupole motion of m Ca and to the extraction of 
transport coefficients from the mean values and dispersions of the quadrupole operator 
Q and its time derivative Q. The results will also be compared to the mean- field limit 
TDHF as well as to the semiclassical limits of VUU and solutions based on the Vlasov 
equation. Furthermore, in Section 4 the specific heat, entropy and free energy of the 
nucleus 40 C a will be evaluated as a function of temperature in the correlated two-body 
theory in comparison to the mean- field limit. A summary and discussion in Section 5 
concludes this study. 



2 Nuclear quantum correlation dynamics (NQCD) 

The dynamical description of the nuclear many-body problem in the non-relativistic 
limit is based on coupled equations of motion for the one-body density matrix p(ll'; t) 
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and the two-body correlation function C2(12, 1'2'; t) in the equal time limit |5^] - [il 
p(ll'; f) = [t(l) + t)]p(ll'; t) - p(ll'; t)[t(l') + C/t(l'; t)} 

+ Tr (2=20 Kl2)c 2 (12, l'2'; t) - c 2 (12, 1'2'; tMl'2')] (1) 

where the two-body density matrix is written as p 2 (12, 1'2'; t) = p 20 (12, 1'2'; t) + 
c 2 (12,l'2';t) with p 20 (12, 1'2'; t) = A 12 p(ll';t)p(22';t) . In (0) t(i) denotes the kinetic 
energy of particle i while i> (ij) is the bare two-body interaction. The time evolution of 
c 2 is determined by 

id/dt c 2 (12, 1'2'; t) = [t(l) + t(2) + C/(l; t) + U(2; t)]c 2 (12, 1'2'; t) 

-08(12, 1'2'; t)[t(l') + t(20 + CT+Cl'; t) + tft(2'; t)] 

+ [gr 2 ^(12)p 2 o(12, 1'2'; f) - p 20 (12, 1'2'; t)v(l'2')Q? 2 ,] (2) 

+ [gr 2 ^(12)c 2 (12, 1'2'; t) - c 2 (12, 1'2'j ^(l^')^] (3) 

+Tr (3= 3 [v(13) AvzAvv ~ v{l'3')AvvA l2 \p{ll'; t)c 2 (23, 2'3'; t) 

+ Tr (3= 30 b(23)^ 23 A' 2 ' - t;(2 , 3 , )^3^i2]p(22 / ; t)c 2 (13, l'3'; t) (4) 

+ Tr (3=3 '){Kl3) + u(23)]c 3 (123, 1'2'3'; t) - c 3 (123, lW; *)Kl'3') + u(2'3')]}. (5) 

In equations (0) - © Trn =i i\ stands for a summation over spin <r(i) and isospin r(i) 
as well as an integration over d 3 ri with i — i', while .Ay is the antisymmetrization 
operator for fermions defined by Aij = 1 — Pij and is the permutation operator 
between particle i and j. In (D) - (|) we have, furthermore, introduced the mean-field 
potential 

U(i; t) = Tr {3=3l) {v(i3)A l3 p(33'; t)} (6) 
and the Pauli-blocking operator for p — p and h — h interactions 

Q= = 1 - 7V (3= 3<) (P,3 + ^3)p(33'; t) (7) 

In order to close ([!]) - @ the three-body correlation term @ will be dropped in the 
following. For an explicit discussion of (|5]) in the context of trace conservation laws we 
refer the reader to Ref. |j4l~| . 



2.1 Time-Dependent Density-Matrix Theory (TDDM) 

Because of technical limitations the set of coupled equations for the one-body density 
matrix p(ll'; t) (0) and the two-body correlation function c 2 (12, 1'2'; t) @ - @) can, at 
present, only be solved accurately for moderate excitations of the nuclear system. In 
this case one can expand both p and c 2 in terms of single-particle states ip a fulfilling 
the TDHF equations, 

(i^-/i(l))tMM) = (8) 
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according to 

p(ll / ;t) = ^n a/3 (t)^(l',t)^(l,t) 
c 2 (12,l'2';t) = CW(*) ^ Q (l,t)^(2,t)C,(l',t)^(2',t), (9) 

oc,P,a',l3' 

where the operator h(i) = t(i) + U(i) in (|8]) is the one-body hamiltonian, i.e. the sum 
of a kinetic energy term and the selfconsistent mean field (|6|) . 

The equations of motion for the occupation matrix n a p(t) and C a p a i^{i) then read 

i ^ n "/3 = XX( aor M^)CW<x - C a6la {^a\v\l35)} (10) 

i -Q^C a p a 'pi = B a p a >pi + P a f3a'f3> + Hapaipi, (11) 



and 



where 

B a (3a>l3> = X! (^1^2MA3A 4 ).4 
A1A2A3A4 

' {(<W - ^aAiX^Aa - ™/3A 2 )™A 3a 'nA 4 /3' ~ n aXi n l3 x 2 (5x 3 a' ~ ™A 3 <*' ) (h 4 /3' ~ ™A 4( 3')} ( 12 ) 

represents the lowest order contribution of collisions in the particle-particle channel 
(Born approximation), while the term P represents the higher order particle-particle 
(and hole- hole) contributions, 

P a /3a'f3' = (^1^2\v\^3^)[^aX 1 ^^X 2 Cx 3 X 4a 'l3' — Sx aa >5x i ^C al 3x 1 X 2 

Ai A2A3A4 

— £aAi"-/3A 2 CA 3 A 4 a'/3' — 5x 2 f3'n a x 1 Cx 4 X 3 f3'a' + ^A 3 a'^A 4 /3'Ca/3AiA2 + ^A 4 /3' n A 3 a' C a f3\ 1 A 2 ] • (13) 

The last term 

Hal3a'/3' = ( Al A2 |l> | A3A4) 

Ai A 2 A 3 A 4 

[^aAi(«A 3 a'C/3A 4 /3'A 2 ~ n A 3 /3'C/3A 4 Ca 3 /3A 2 /3' — '"-A 4 /3'C , A 3 /3a'A 2 ) 

+^/3A 2 (^A 4 /3'CQA 3 a'Ai — ™A 4 a'CaA 3 /3'Ai ~ ^A 3 /3'C , QA 4 a'A 1 ~ '^A 3 a'CaA 4 A 1 /3') 

— ^/3'A 4 ( ra /3A 2 CaA 3 a'A 1 ~ ra aA 2 C/3A 3 a'Ai ~ ^/3Aj C«A 3 a'A 2 — n a\ 1 Cf3\ 3 \ 2 a) 

— 5q'A 3 (^qAiC/3A 4 /3'A 2 — "-aA 2 CaA 4 /3'A 2 ~ ^aA 2 C/3A 4 /3'Ai — ^/3A 2 CaA 4 Ai/3')] (14) 

is the contribution to the equations of motion in the particle-hole channel. The set of 
coupled equations (|I0| ) and ( |TT| ) provides a non-perturbative description of large am- 
plitude nuclear motion denoted as Time-Dependent Density-Matrix (TDDM) theory 
which takes into account collisions among nucleons to all orders in the interaction and 
thus superseeds the quantum approaches in Refs. [^, |42| that are limited to the Born 
approximation, i.e. retaining only B a i vha p a ipi in Eq. fllT]). We note that the theory 
fulfills the conservation of particle number, momentum, and energy |34| . 
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2.2 Dispersion of one-body operators 

Since TDDM is a consistent two-body theory it allows to evaluate one-body as well 
as two-body observables. Of specific interest are the mean values and the fluctuation 
properties of collective one-body operators, for which we give the microscopic expres- 
sions. Note that contrary to Refs. [E3l [37], B3[ no explicit equations of motion for a 



collective degree of freedom have to be set up since Eqs. (10) and (11) already describe 
the full dynamics. 

The quadrupole operator (as well as any one-body operator O) is defined as 

O = J2(oc\0\(3)aia p , (15) 

where a+ and ap are creation and annihilation operators of a nucleon. Our definition 
of the quadrupole operator Q2(r) = Q(r) is explicitly 

Q(r) = l^J(2z 2 -a; 2 -y 2 ) (16) 

which gives for the time derivative 

— i / 5 1 

<5( r > P) = ~^[Q,h] = W -^7^( 2z Pz + 2 P* Z ~ x Px ~ PxX - yp y - VyV)- (17) 

The inertia mass operator of the quadrupole motion, M(r), is 

= -^K?00, [Qto, h\\ = -l[g(r), Q(r, p)] = JL ( 4 , 2 + x 2 + y 2 ) . (18) 

Eq. (|18D presents not only an expression for the quadrupole mass operator but also a 
more general commutation relation for the quadrupole operator Q(r) and its velocity 
operator Q(r,p), in which the constant ih is replaced by the operator ih/M(r). This 
is equivalent to a collective momentum operator P = M(r)Q(r, p). 

From the commutation relation (|IBD one can derive the corresopnding uncertainty 
relation as, 

^^(I<mW>) 2 ' (19) 

or 

°lMy Q >(^\ (20) 

with the definition of the quadrupole mass 1/Mq = (1/M(r)), and with Oq and cr 2 ^ 
defined by Eq. (p2|) below. Eq. (|20|) will be checked as a boundary condition in our 
following numerical calculations explicitly. 

The mean value and the dispersion of a one-body operator O are given by: 

(O) = J2W\0\P){atafs) = J2( a \°\P) n 0<* (21) 
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°l= E HO|/3)(e/|O|/30(a^a>^)-(O) 2 

a,a',/3,/3' 

a,a',/3,l3' 

= Y.^\0\f3){f3\0\l3')np la + £ (a|0|/3)(a'|0|)9 , >(C / j /9 / aa /-n j8a /n / y a ) 

a,/3,/3' a,a',/3,f3' 

= E(«|C 2 |/3>^+ E (alOI^^'IOI^^w-^m^), (22) 

where the matrix elements and Co^/^ depend on time explicitly. In case of 
TDHF we simply have Cpp aa i = and n ai g = b a pn a where n Q denotes the occupation 
probability of the single-particle state ip a . This gives 

(O) =J2(a\0\a)n a , (23) 

a 

a 2 = J2(^\0\(3)(l-n p )(f3\0\a)n a 

a,/3 

= J2(a\0 2 \a)n a - ^2(a\O\/3)(j3\O\a)n0n a . (24) 

We note that the summation over (3 in the first terms of Eqs. ([22] ) and (^4j) should be 
taken in the whole Hilbert space (J2p \P >< P\ — 1) so that it yields an identity in the 
whole Hilbert space. This observation is important since the numerical calculations 
are carried out in a truncated Hilbert space. 



2.3 Technical approximations 

In our actual calculations we approximate the interaction v appearing in the mean-field 
potential (|6]) by a Bonche-Koonin-Negele (BKN) |P| force and use v (12)=Vo<5(ri — r 2 ) 
for the residual interaction appearing in (|I1^ -(|T4"D with Vo = —300 MeV fm 3 for tech- 
nical reasons (see below) as in Refs. J7], ^ . We have solved the coupled equations 
(|10"D and ([□]) for the case of the isoscalar monopole or quadrupole vibration of 40 Ca as 
a function of temperature since the dipole mode was found in Ref. [7] to be dominated 
by mean-field dynamics. To carry out calculations at various temperatures we describe 
the initial occupation numbers by appropriate Fermi distributions (characterized by a 
temperature T) and propagate the system in time. The set of single-particle levels used 
in the calculations include the Is, lp, 2s, Id, 2p, and If shell. The correlated state, 
which is generated by an adiabatic switching on of the residual interaction j|5, 



is then boosted at a certain time (typically after 2 • 10~ 22 s) by applying appropriate 
phase factors to the wavefunctions ipp proportional to a strength factor a. In this way 
the system acquires a well defined collective energy proportional to a 2 . The set of 
equations (|S[) - ( |TTD is then integrated in time with timestep-size At = 0.5 • 10 _23 s. We 
note that the total energy and number of particles are conserved well throughout the 
timescales of interest here. 
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2.4 Specific heat, entropy, and free energy 

The energy of the system reads for the present Skyrme interaction as well as the residual 
S— interaction: 

E = J e(r) dh = J d 3 r [|-(r n + r p ) + ^(2p nPp + i(p£ + p 2 p )) 

+ j(PnPp + Pnp 2 p )} + ^-[E y (p n , p n ) + E y (p p , Pp )} + V V E y (p ni Pp ) 

+ ?- [ [ dhdh> + ^ I C 2 (r, r; r, r) dh, (25) 



r — r 



where 



E y ( Pq , Pq ') = I I dhdh' 6XP( Hr ) , r/|) ft(rW(0 (26) 



p | r — r 



results from the finite range Yukawa force The total energy of the system - after 
adiabatically switching on the residual interactions - depends on the initial thermal 
excitation energy characterized by a temperature T, from which one can calculate 
numerically the specific heat as 

cm = (27) 

The entropy for the many-body system, which should not be mixed up with the one- 
body entropy, then is given by 

S(T) = J q T ^pdT' (28) 

and the free energy follows as 

F{T) =E-TS. (29) 
Furthermore, the fluctuation in energy of the system can be evaluated as, 

a 2 H = ({H - (H)) 2 ) = (H 2 ) — E 2 = T 2 C V , (30) 

which can be considered as an order parameter for a phase transition. 



2.5 The semiclassical limit 

As shown in detail in Ref. [|J the coupled equations of motion (JJ) - @, i.e. excluding 
(4), transform to a Vlasov-Uehling-Uhlenbeck (VUU) equation in the semiclassical 
limit when adopting an on-shell (Markov) approximation for the iViV-interactions and 
expanding up to the first order in the derivative operator fr -1 V r • V p with respect to 
the particle phase-space variables r and p. This gives the familiar VUU equation (for 
momentum-independent mean-fields) 

{^7 + ^ ■ V r - V r C/(r) • V p }/(r, p; t) = I coU (r, p; t), (31) 
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where the collision term is given by Q 



4 f f dcr 
Icoii(r,p;t) = j— 5 J d 3 P2 J dVL K 2 |— (p,p 2 ) 

x (/(r, p 3 ; t)f(r, p 4 ; t)(l - /(r, p; t))(l - /(r, p 2 ; t)) 
- /(r, p; t)/(r, p 2 ; t)(l - /(r, p 3 ; *))(1 - /(r, p 4 ; t)). (32) 
In Eq. ( |3~2"D dcr/dfi denotes the differential cross section for the on-shell scattering 



process 1 + 2 — > 3 + 4 while |t>i 2 | is the relative velocity of the colliding nucleons; 
the outgoing momenta p 3 and p 4 are fixed by energy- and momentum conservation 
except for the scattering angle Q in their center-of-mass. When neglecting the collision 
term (0), the Vlasov equation (|3lD is equivalent to the statement, that the total 
time derivative of the one-body phase-space density /(r, p;£) vanishes, i.e d/dtf = 
0. The Vlasov equation thus corresponds to the semiclassical limit of TDHF. The 
VUU equation ( [311) is conventionally || |[] solved within the testparticle method and 
has been exploited extensively for the nonequilibrium description of proton-nucleus 
or nucleus-nucleus collisions at intermediate H and high energies ^7|. Here we will 



compare the results of the semiclassical VUU or Vlasov equation with those for the 
quantum mechanical theories TDDM and TDHF, respectively. 



2.6 The damped harmonic oscillator in quantum physics 

For later comparison and interpretation of the numerical results we recall the equa- 
tions of motion for a damped harmonic quantum oscillator (QHO) with respect to 
fluctuations in the collective coordinates Q and Q (cf. Refs. [^, [3^ 



Jfl = 2a QQ ( 33 ) 

2 ^QQ + TT°Q (34) 



d 2 . o 2 ^ 



L Q 

where the fluctuations are defined as 



= -2 lQ *< + — a QQ + Dqq(T), (35) 



° 2 Q = (Q 2 ) - 


(Q) 2 ; 


qq = (QQ) - 


(QHQ) 


°l = (Q 2 ) - 


(Q) 2 - 



(36) 

In Eq. ( p4|) Q denotes the collective frequency, 7q a collective friction coefficient, Mq 
the collective mass parameter, whereas Dqq(T) is the diffusion coefficient given by 

DQQ{T) ^^ mth m (37) 
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Note that only for large temperatures T 3> hQ the Einstein relation 

D QQ = 2 lQ T/M Q (38) 

is obtained. In view of (small amplitude) nuclear collective quadrupole motion we have 
T < Ml such that relation (|3~8"D does not apply. When using P = MqQ instead of Q 
the diffusion coefficient ( fTFD or (|38D picks up an additional factor Mq. 

3 Numerical results 

As a first test we study the 'hole-strength' distribution of 40 C a within the nonperturba- 
tive limit TDDM by integrating the coupled equations (|), (]TD|) , and ([□]) in time with 
the initial condition C a p a ipi = 0, which corresponds to the respective Hartree-Fock 
groundstate. The correlated state is then generated in line with the Gell-Mann-Low 
theorem by adiabatically switching on the residual interaction as 

V(t) = V [l-exp(-t/r)], (39) 

where r ~ 5- 10~ 22 sec was found to be sufficiently large to achieve convergence (cf. also 
Ref. P§|). The numerical solutions for the occupation numbers n a (t) deviate from 
and 1 and slightly oscillate in time. By averaging over time (for t > 10~ 21 sec) we obtain 
average occupation numbers (n a ) which are shown in Fig. 1 as a function of the time- 
averaged single-particle energies e a . The solid line presents a smooth fit through the 
discrete points for orientation. The latter one is remarkably close to the experimental 
data from Ref. f|8| which indicates that the proper strength of the residual correlations 
is well described by TDDM for Vq = —300 MeV fm 3 , which we consider as fixed from 
now on. 

We note in passing that the average occupation numbers have no direct counter 
part in the semiclassical limit since the s.p. energy does not enter the VUU equation 
explicitly, but only its space- and momentum derivatives. In analogy to the mean-field 
limit the nucleon phase-space distribution at T = is given by a local Thomas-Fermi 
distribution (for fixed spin and isospin) 

/o(r J p;t) = e(p i r(r)-|p|) (40) 

with 

p F (r) = (^p(r)) 1/3 (41) 
where p(r) denotes the nucleon density distribution. 

3.1 Damping of giant quadrupole motion 

Whereas in mean-field theory the frequency of collective modes can be rather well 
determined, its damping width appears to be dominated by higher order NN interac- 
tions. As was found in Ref. the role of 2p — 2h matrix elements is very pronounced 
for giant quadrupole motion, which we analyse here in more detail. 
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The time evolution of the quadrupole moment Q(t) and that of Q(t) - which carries 
a similar information - are presented in Figs. 2 and 3 for TDHF (upper parts) and 
TDDM (lower parts) for an isoscalar quadrupole excitation of 20 MeV of 40 Ca at an 
initial temperature of T=l, 3, and 5 MeV. 

Whereas there is only small damping in the case of TDHF, the collective motion is 
strongly damped in the limit TDDM demonstrating the decisive role of residual NN 
interactions. 

The questions arises, to what extent the damping of the quadrupole motion can be 
attributed to on-shell particle-particle collisions as inherent in the semiclassical VUU 
limit. To this aim we show in Fig. 4 the time evolution of the quadrupole moment in 
momentum space, 

Q * {t) = (2W Mp M -rf -rf)/( r .p;*). ( 42 ) 



where /(r, p;i) results from the VUU equation ([31]) in either the mean-field limit 
(denoted by Vlasov) or the VUU limit including the on-shell two-body collisions (VUU). 
As in case of the quantal approaches there is only a minor damping in the Vlasov case 
whereas a rapid decay of Q^if) is seen in the VUU limit. Unfortunately, no precise 
temperature T can be attributed to the semiclassical calculations due to numerical 
reasons since the testparticle method involves numerical fluctuations in the collective 
coordinate as well as in the total energy. From the size of the numerical fluctuations 
we estimate the 'excitation' temperature to be in the order of 3-5 MeV. 

Since the VUU approach includes only on-shell particle-particle collisions and leaves 
out particle-hole interactions, which are important for the damping of shape vibra- 
tions |JT|, it yields a moderate damping in between the limits TDHF and TDDM. 



From the height and position of the maxima (and minima) we can extract the average 
collective frequency Q in all cases as well as an average friction constant Tq using 
Q(t) ~ exp(— Tqt/2h). The latter quantity is shown in Fig. 5 for initial nuclear tem- 
peratures T of 0, 1, 2, 3, 4, 5, and 6 MeV for TDHF and TDDM where the error bars 
indicate the uncertainty of extraction from the time-dependent signal Q(t) or Q(t). 

The actual values are found to be approximately independent on temperature for 
TDHF, whereas the damping drops with T for TDDM. Fig. 5, furthermore, shows 
that the damping achieved in TDHF compares well with the semiclassical Vlasov limit 
for T ~ 3 — 5 MeV and the additional damping introduced by collisions in VUU only 
provides about half of the friction obtained from TDDM at low temperature (T « 0). 
This result indicates the essential role played by the particle-hole correlations (long- 
range correlations) in the damping of shape and surface vibrations. As discussed in 
IS HH a sizeable damping arizes from the particle- hole (p-h) interaction matrix H a p a /pt 
014]) that is not considered in the limit with nonperturbative p-p and h-h scattering 



(e.g. VUU like models or in Ref. fl33[). 

The actual mag nitude of T Q « 5 - 6 • lO 21 ^ 1 at T ^0 from TDDM compares well 
with the friction constant extracted experimentally from the width of the GQR for 
m Ca and indicates the importance of p-h residual matrix elements. 

The dropping of the width Tq with temperature for TDDM in case of isoscalar 
quadrupole motion can be attributed to the fact that the long-range correlations de- 
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crease in magnitude with temperature. In order to demonstrate this dependence we 
show in Fig. 6 the dependence of the correlation energy on temperature T which in 
magnitude decreases by a factor of 2 from T = to T = 6 MeV. Thus the 'long-range 
correlations melt with temperature' and the width of the giant quadrupole resonances 
should decrease with T. This is opposite to the role of on-shell 2-body collisions which 
are blocked in the VUU limit at T=0 and increase with phase space as ~ T 2 . Off-shell 
collisions allow for a faster increase of phase-space with temperature as demonstrated 
in Ref. [33"||, however, the functional dependence of the collisional width depends cru- 



cially on the multipolarity of the mode. Thus dominantly the particle-hole correlations 
induced by (|l4]) generate sizeable density fluctuations at low temperature T which lead 
to a coupling of different collective modes and a redistribution of their strength. 

3.2 Dispersion of giant quadrupole motion 

In the quantum theories the dispersion of one-body operators can be evaluated explic- 
itly and the effect of 2-body correlations on these 2-body operators can be analyzed in 
detail. We note that a collective mass parameter can be extracted from the calculations 
via 

E„u - ^(Q)L X - ^ 2 (Q) 2 max (43) 

and leads to Mq w 2.0MeV/(fm 2 c 2 ) which is practically equivalent with the value 
of the direct calculation from Eq. (|18|). The calculated dispersions Oq(£) and o%{t) 
according to ( f22"| ) are found to oscillate in phase with the quadrupole moment (Q)(t) 
or (Q)(t) (cf. Fig. 7). When using the same initial conditions cr^(t = 0) and Oq(t = 0) 
as well as the collective frequency Q and 7q = Tq from the TDDM calculations, the 
quantum harmonic oscillator (QHO, Eqs. (33)-(35)) gives comparable solutions for 
Oq(£) and cr 2 ^{t) which, however, are damped out faster in time (dotted lines in Fig. 
7). In case of TDDM or TDHF the solutions still oscillate around some average value. 

The asymptotic values of <Tq and a? (for t — > oo) are obtained by averaging over 
the final time interval (typically for t > 7.5 ■ 10~ 22 sec) and shown in Fig. 8 as functions 
of the initial temperature T. Whereas o 2 q and a 2 ^ only slightly varies with temperature 
for the QHO, the variation with T is more pronounced for TDDM (solid lines) and 
especially TDHF (dashed lines), which reflect the intrinsic dynamics. Thus model 
studies based on the QHO equations of motion should be taken with care especially at 
high temperature |43[ 



The product MqU^ctq according to (|20|) must obey the uncertainty relation in all 
quantum mechanical approaches. This is demonstrated in Fig. 9 for TDDM (solid 
line), TDHF (dashed line) and the QHO (dotted line) as a function of temperature T. 
Note, that h 2 « 0.4 (MeV 10" 21 s) 2 . 

Since \Mqo 2 ^ together with the potential energy \MqV?Oq should reflect the 



equipartition theorem in case of damped harmonic motion, [23 



E Q (T) = \M Q a 2 Q + l -M Q Va 2 Q = E Q (0) + T, (44) 
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and thus be a linear function of the temperature T according to Eq. (PD, we have 
plotted 

AE Q (T) = E Q {T) - E Q (0) (45) 

in Fig. 10 for all three cases. Here only the full two-body theory TDDM follows 
equipartition for T > 1 MeV, i.e. AEq(T) = T, whereas TDHF as well as the QHO 
significantly overestimate or underestimate the dash-dotted line, respectively. The 
off-set of ~ 1 MeV for TDDM at T = so far is not understood. 



4 Specific heat, entropy and free energy 

Apart from the crucial role of long-range correlations in the damping of the collective 
quadrupole motion demonstrated in Section 3, the two-body correlations are expected 
to modify also the thermodynamic properties relative to a mean-field approach. To 
investigate this question we calculate the total energy of our test-nucleus 40 Ca (without 
any collective excitation) as a function of the initial temperature from T— to T = 6 
MeV. We note that due to two-body correlations the binding energy is lowered by pa 
20 MeV (i.e. -382 MeV) for TDDM as compared to TDHF (- 362 MeV). 

In Fig. 11 we show the total energy for TDDM (full squares; shifted by 20 MeV) 
and TDHF (full circles) as a function of temperature (upper part) as well as the specific 
heat ( p?D (lower part). The total energy increases slightly faster with T for TDDM as 
for TDHF due to a decrease in the correlation energy with T (cf. Fig. 6). Whereas 
the increase of the energy is monotonous for TDHF, there is a more rapid variation of 
the total energy for TDDM at T pa 4 MeV. This is seen more clearly in the specific 
heat (lower part), where Cy shows a sudden bump around T « 4 MeV, whereas the 
specific heat is flat as a function of T for TDHF. 

Since conventionally maxima in the specific heat are attributed to phase transitions, 
where the fluctuations in energy become very large, we have evaluated the quantity 
< H 2 > — < H > 2 according to fl30l) . As seen from Fig. 12 the fluctuations in the 
total energy also show a bump at T pa 4 MeV for TDDM (solid triangles) whereas they 
stay rather flat for TDHF (open circles) in this temperature region. 

The origin of this bump is not due to a rapid change of the correlation energy (cf. 
Fig. 6) but reflects the melting of single-particle effects in case of TDDM. This is 
demonstrated in Fig. 13 where we show the single-particle density p{x = y = 0; z) as 
a function of the coordinate z for temperatures from T=0 to T = 6 MeV (upper part) 
as well as the single-particle potential (lower part). The pronounced maxima (for z ~ 
0) in case of low temperature T < 3 MeV no longer appear for T > 4 MeV contrary 
to the case of TDHF. We attribute this behaviour to the melting of shell structure by 
the two-body correlations in case of TDDM which is accompanied by a larger specific 
heat as well as larger fluctuations in the total energy (cf. Figs. 11 and 12). 

We finally note that the role of two-body correlations also shows up in the entropy S 



fl2~8|) and the free energy F (p9|) of the nucleus. The latter two quantities are displayed 
in Fig. 14 as a function of temperature for TDDM and TDHF. Whereas the two-body 
correlations lead to a decrease of the free energy F relative to the mean-field result, 
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they enhance the entropy of the system due to a stronger dynamical mixing of single- 
particle configurations (Slater-determinants). However, the relative changes turn out 
to be quite moderate except for a shift of the free energy by ~ 20 MeV for TDDM due 
to the correlation energy. 



5 Summary and discussion 

In this work we have presented microscopic calculations for the nucleus m Ca within 
a nonperturbative two-body correlation approach which superseedes TDHF as well as 
the 'traditional' VUU model || |J by including explicit p-h interactions as well as 
off-shell collision processes to all orders. Note that more recent extensions of mean- 
field theory |33|, |42| include off-shell collisions, however, are still restricted to the Born 
approximation. 

By studying the decay of collective quadrupole excitations we have extracted trans- 
port coefficients for friction and diffusion and analyzed the fluctuation properties of 
the theory which fulfills the equipartion theorem for the fluctuations in the collective 
velocity (for T > 1 MeV) contrary to the mean-field limit TDHF. This study has also 
shown the essential role played by the particle-hole (long-range) correlations in the 
damping of shape and surface vibrations in nuclei at lower temperature whereas for 
T > 3-4 MeV the particle-particle collisions take over and the system might be conve- 
niently described within the semiclassical VUU limit. The latter phenomena is due to 
a decrease of the correlation energy with temperature and an increase with phase-space 
of on-shell two-body collisions with T. 

The temperature dependence of the entropy S and the free energy F is found not 
to be modified substantially due to two-body correlations relative to the mean-field 
limit. However, a small bump in the specific heat Cy and in the fluctuations in energy 
shows up for T « 4 MeV, which does not appear in mean-field theory. We attribute 
this effect to a melting of shell structure due to two-body correlations. 

In view of the detailed microscopic analysis we expect that nucleus-nucleus collisions 
above 20 A MeV should be adequately described in the semiclassical VUU limit as far 
as stopping power and single-particle observables are concerned. Multifragmentation 



phenomena, however, might require the dynamics of long-range correlations p6| , [49 
especially at low nuclear density. Mean-field dynamics here provide a glance at the 
complexity of the shape evolution KO, BTJ but should not be adequate to describe 



nuclear condensation phenomena or the growth of instabilities according to the analysis 
in this work. First steps in exploiting Eqs. (1) - (4) for the description of nucleus- 
nucleus collisions have been made in Ref. 115 211 . 



References 

[1] A. A. Abrikosov, L. P. Gorkov and I. E. Dzyaloshinski, Quantum field theoretical 
methods in statistical physics, (Pergamon, Oxford, 1965) 



13 



[2] A. B. Migdal, Theory of finite Fermi systems and application to atomic nuclei, 
(Wiley Interscience, New York, 1967) 

[3] G. F. Bertsch and S. Das Gupta, Phys. Rep. 160 (1988) 189 

[4] W. Cassing, V. Metag, U. Mosel, and K. Niita, Phys. Rep. 188 (1990) 363 

[5] A. Bonasera, F. Gulminelli, and J. Molitoris, Phys. Rep. 243 (1994) 1 

[6] J. Wambach, Rep. Prog. Phys. 51 (1988) 989; S. Drozdz, S. Nishizaki, J. Speth, 
and J. Wambach, Phys. Rep. 197 (1990) 1 

[7] F.V. De Blasio, W. Cassing, P. F. Bortignon, R. A. Broglia, and M. Tohyama, 
Phys. Rev. Lett. 68 (1992) 1663 

[8] S. Ayik and C. Gregoire, Phys. Lett. B212 (1988) 269; Nucl. Phys. A513 (1990) 
187 

[9] J. Randrup and B. Remaud, Nucl. Phys. A514 (1990) 339 

[10] G. F. Burgio, Ph. Chomaz and J. Randrup, Nucl. Phys. A529 (1991) 157; Phys. 
Rev. Lett. 69 (1992) 885 

[11] J. Randrup, Nucl. Phys. A545 (1992) 47c 

[12] S. Ayik, E. Suraud, M. Belkacem, and D. Boilley, Nucl. Phys. A545 (1992) 59c 

[13] P. G. Reinhard and E. Suraud, Nucl. Phys. A545 (1992) 59c; Ann. of Phys. 216 
(1992) 98 

[14] A. Bonasera, F. Gulminelli, and P. Schuck, Nucl. Phys. A545 (1992) 71c 

[15] Ph. Chomaz, G.F. Burgio, and J. Randrup, Nucl. Phys. A540 (1992) 227 

[16] Y. Abe, S. Ayik, P. G. Reinhard, and E. Suraud, Phys. Rep. 275 (1996) 49 

[17] M. Germain et al., Phys. Lett. B437 (1998) 19 

[18] T. Maruyama, A. Onishi and H. Horiuchi, Phys. Rev. C45 (1992) 2355 

[19] A. Bracco et al., Nucl. Phys. A519 (1990) 47c 

[20] G. Enders et al., Phys. Rev. Lett. 69 (1992) 249 

[21] J. J. Gaardhoje, Annu. Rev. Nucl. Part. Sci. 42 (1992) 483 

[22] E. Ramakrishan et al., Phys. Rev. Lett. 76 (1996) 2025; Phys. Lett. B383 (1996) 
252 

[23] W. Cassing and W. Ndrenberg, Nucl. Phys. A401 (1983) 467 

[24] P. F. Bortignon, A. Bracco, D. M. Brink, and R. A. Broglia, Phys. Rev. Lett. 67 
(1991) 3360 

14 



[25] Ph. Chomas, M. Di Toro and A. Smerzi, Nucl. Phys. A563 (1993) 509 
[26] A. Smerzi, M. Di Toro and D. M. Brink, Phys. Lett. B320 (1994) 216 
[27] C. Toepffer and P.-G. Reinhard, Ann. Phys. 181 (1988) 1 

[28] K. Gutter, P.-G. Reinhard, K. Wagner, and C. Toepffer, Ann. Phys. 225 (1993) 
339 



[29] M. Di Toro, V. M. Kolomietz, and A. B. Larionov, |nucl-th/9807070| 



[30] P. F. Bortignon, M. Braguti, D. M. Brink et al, Nucl. Phys. A583 (1995) 101c 

[31] G.F. Bertsch, PF. Bortignon, and R.A. Broglia, Rev. Mod. Phys. 55 (1983) 287 

[32] N. D. Dang, K. Tanabe and A. Arima, RIKEN-AF-NP-285, 1998; RIKEN-AF- 
NP-296, 1998 

[33] D. Lacroix, Ph. Chomaz and S. Ayik, Phys. Rev. C58 (1998) 2154 

[34] M. Tohyama and M. Gong, Z. Phys. A332 (1989) 269; M. Gong and M. Tohyama, 
Z. Phys. A335 (1990) 153 

[35] M. Gong, M. Tohyama and J. Randrup, Z. Phys. A335 (1990) 153 

[36] W. Cassing, A. Peter and A. Pfitzner, Nucl. Phys. A561 (1993)133 

[37] S. Ayik and J. Randrup, Phys. Rev. C50 (1994) 2947 

[38] S. J. Wang and W. Cassing, Ann. Phys. 159 (1985) 328 

[39] W. Cassing and S. J. Wang, Z. Phys. A328 (1987) 423; Z. Phys. A337 (1990) 1 

[40] W. Cassing and A. Pfitzner, Z. Phys. A337 (1990) 175 

[41] W. Cassing and A. Pfitzner, Z. Phys. A342 (1992) 161 

[42] D. Lacroix, Ph. Chomaz and S. Ayik, nucl-th/ '9902050 

[43] W. Wen, P. Chau Huu-Tai, D. Lacroix, Ph. Chomaz, and S. Ayik, Nucl. Phys. 
A637 (1998) 15 

[44] P. Bonche, S. E. Koonin, and J. W. Negele, Phys. Rev. C13 (1976) 1226 

[45] A. Pfitzner, W. Cassing, and A. Peter, Nucl. Phys. A577 (1994) 753 

[46] A. Peter, W. Cassing, J. M. Hauser and A. Pfitzner, Nucl. Phys. A573 (1994) 93 

[47] W. Cassing and E. L. Bratkovskaya, Phys. Rep. 308 (1999) 65 

[48] L. Lapikas, Nucl. Phys. A553 (1993) 297c 

[49] H. Heiselberg, C. J. Pethick, and D. G. Ravenhall, Nucl. Phys. A519 (1990) 279c 

15 



[50] C. Jung, W. Cassing, U. Mosel, and R. Y. Cusson, Nucl. Phys. A477 (1988) 256 

[51] D. Lacroix and Ph. Chomaz, Nucl. Phys. A641 (1998) 107; Phys. Rev. C58 (1998) 
1604 

[52] Jian-Ye Liu, Shun- Jin Wang et al, Nucl. Phys. A604 (1996) 341 



16 




i 1 1 1 1 1 r 

-15 -10 -5 



s.p. energy (MeV) 



Fig. 1: The time averaged occupation numbers n a for 40 Ca as a function of the s.p. energies 
e a at T = MeV in the limit TDDM. 
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Fig. 2: Time dependence of the quadrupole moment Q2{t) = Q(t) for an isoscalar excitation 
of 40 Ca of 20 MeV in the limits TDHF (upper part) and TDDM (lower part) at initial 
temperatures of 1,3, and 5 MeV. 
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Fig. 3: Time dependence of the quadrupole velocity Q2(t) = Q(t) for an isoscalar excitation 
of 40 Ca of 20 MeV in the limits TDHF (upper part) and TDDM (lower part) at initial 
temperatures of 1, 3, and 5 MeV. 
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Fig. 4: Time evolution of the quadrupole moment in momentum space (42) in the limits VUU 
(solid line) and Vlasov (dashed line) 
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Fig. 5: Extrapolated values for the friction coefficient Tq as a function of the initial temper- 
ature T in the limits TDDM, TDHF, VUU, and Vlasov. 
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Fig. 7: The fluctuations ag(t) (upper part) and (lower part) for 40 Ca excited collectively 

by 20 MeV in the limit TDDM (full lines) and TDHF (dashed line) in comparison to the 
results from the damped quantum harmonic oscillator (QHO; dotted lines). 
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Fig. 8: Asymptotic fluctuations in the quadrupole degree of freedom Oq(oo) (upper part) and 
in the quadrupole velocity a^(oo) (lower part) in the limits TDDM, TDHF and the quantum 
harmonic oscillator QHO. 
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Fig. 9: Test of the uncertainty relation for TDDM (solid line), TDHF (dashed line) and the 
QHO (dotted line) function of T. 
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Fig. 10: Test of the equipartition relation for AEq(T) (45) (straight dash-dotted line) for 
the collective fluctuations as a function of temperature in the limits TDDM, TDHF and the 
quantum harmonic oscillator QHO. 
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Fig. 11: The total energy (upper part) and the specific heat (27) (lower part) of the system 
40 Ca as a function of temperature in the limits of TDHF (open circles) and TDDM (full 
squares) . 
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Fig. 12: The fluctuation in energy (30) of the system m Ca clS Si function of temperature in 
the limits of TDHF (open circles) and TDDM (full triangles). 
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Fig. 13: The nuclear density (upper part) and the single-particle potential (lower part) 
for 40 Ca as a function of the coordinate z for temperatures T— 1,2,3,4,5,6 MeV in the 
limit TDDM. 
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Fig. 14: The entropy S (28) (upper part) and free energy F (29) (lower part) for 40 Ca as a 
function of temperature for TDDM and TDHF. 
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